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We study possible restoration patterns of chiral symmetry in a generalized hidden local symmetry 
model, which is a low energy effective theory of QCD including pseudo-scalar, vector and axial-vector 
mesons. We derive Wilsonian renormalization group equations and analyze the running couplings 
and their fixed points at the chiral restoration point. We find three types of the chiral restoration, 
which are classified as the standard, vector manifestation and intermediate scenarios, respectively. 

It turns out that the rho and A\ meson become massless and their decay into pion is suppressed 
in all the restoration patterns. The each restoration scenario violates or fulfills the vector meson 
dominance at the critical point in a different manner, which may reflect on the contributions from 
the pion to the dilepton spectrum. 
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I. INTRODUCTION 

Dynamical chiral symmetry breaking is the most remarkable nonperturbative phenomenon in QCD together with 
confinement. One of the most important goals in hot and/or dense QCD is to understand how chiral symmetry 
restores at high temperature and/or density, which is commonly believed to be true. Then, it is crucially important 
to understand how hadrons change their properties in hot and/or dense matter, since it is what one observes in 
the actual experiments. Along these lines there have already been many theoretical as well as experimental works 

iiiiiiii 

Theoretically, the properties of pseudo-scalar mesons are determined by the chiral symmetry through low-energy 
theorems since pseudo-scalar mesons are considered to be approximate Goldstone bosons, while the properties of other 
mesons are not directly controlled by the symmetry alone. An interesting possibility, the vector manifestation (VM), 
has recently been proposed by Harada and Yamawaki ;8], in which the vector meson becomes the chiral partner of 
the pseudo-scalar meson when the chiral symmetry is restored and provides a theoretical basis for Brown-Rho scaling 
{§j, i.e. the dropping of the vector-meson mass in hot and/or dense matter. The VM is based on the hidden local 
symmetry (HLS) by Bando et al. PI IH1313 0, which is a symmetry not manifest in the QCD lagrangian but 
is assumed to be dynamically generated, and vector mesons are introduced as gauge bosons associated with the HLS. 
Though there had been a long history before the work of Bando et al. [Hj UJJ LLa- UJi |2(j , in trying to incorporate 
vector mesons as gauge bosons in the low-energy effective field theory, the HLS made it possible to include vector 
mesons in a systematic loop expansion of the chiral perturbation theory )22H . 

In the broken phase of the chiral symmetry the pseudo-scalar and axial-vector mesons are mixed with each other 
while scalar and vector mesons are pure, when mesons are classified according to the representation associated with 
light-front chiral charges [271 l28t l2fl 13011 . This mixing is expected to vanish in the symmetric phase and the mesons 
form chiral doublets. It is usually believed that the scalar meson becomes the chiral partner of the pseudo-scalar 
meson. This is called the ‘standard scenario’. There is, however, another possibility that the longitudinal part of the 
vector meson becomes the chiral partner of the pseudo-scalar meson. This is the ‘VM scenario’ mentioned above. 
Therefore, it is an interesting question, if the VM scenario is realized in the real world or not. The VM is claimed to be 
realized in the HLS model at large number of flavors in ref, fa l23l 123) by Harada and Yamawaki, at high temperature 
in ref. |2^| by Harada and Sasaki and at high density by Harada, Kim and Rho in ref.[26|] . There argument goes 
as follows. The parameters of the HLS are determined from QCD by the Wilsonian matching proposed in ref. 1241 . 
namely by matching the vector and axial-vector current correlators in the HLS model with those in the operator 
product expansion (OPE) of QCD at the matching scale, A. Then, they obtain the parameters of the HLS at lower 
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physical energy by letting the parameters run from the scale, A, according to the renormalization group equation 
(RGE). Now, when the chiral symmetry is restored, the vector and axial-vector current correlators must coincide with 
each other. When the parameters of the HLS run according to the renormalization group equation, this condition 
is always satisfied while the pion decay constant, which is non-zero at the matching scale, A, goes to zero at lower 
energy, i.e. on-shell. In the HLS model, the condition that the vector and axial-vector current correlators coincide 
together with the vanishingpion decay constant leads to the VM, because only the pseudo-scalar and vector mesons 
are included in the model p]. Thus, it was concluded that the VM is realized in the symmetric phase. 

These works, however, are not fully satisfactory in a sense that the model of the HLS includes only the pseudo-scalar 
and vector mesons but neither scalar nor axial-vector mesons. Namely, from the very beginning the possibility of 
realizing the standard scenario seems to be excluded. Thus, a natural question is what happens if one includes the 
axial-vector and scalar mesons into consideration. The purpose of this paper is to give a partial answer to the question 
by employing the generalized hidden local symmetry (GHLS) model 11. l3ll|. which includes the axial-vector meson 
in addition to the pseudo-scalar and vector mesons. We analyze the renormalization group equation in the GHLS 
model and investigate which scenario is realized in the chiral symmetric phase. It should be noted that importance of 
the axial-vector meson is stressed from the viewpoint of the large n — Ai mixing and also of recent results of STAR 
collaboration phi 33|. 

The paper is organized as follows. The concept of GHLS is introduced in sec. nn The Wilsonian matching for the 
GHLS model is performed in sec. m The renormalization group equations and their fixed points are given in sec. 
m Possible restoration patterns are examined in sec. El Sec. lVU is devoted to the summary of the paper. The actual 
calculation to obtain the coefficients of the renormalization group equations is extremely complicated and only the 
results are shown in sec. HY1 We provide the minimum necessary information in order to reproduce the results in 
Appendix, still which is very long. 


II. HIDDEN LOCAL SYMMETRY 

The concept of the hidden local symmetry is that a symmetry, which is not manifest in the fundamental theory, is 
dynamically generated in the effective theory. The QCD with Nf massless quarks and its low-energy effective theory, 
the non-linear sigma model shares the same global symmetry, Ggiobai = SU(V/)l x SU(AT/)r, under which the pion 
field U is transformed as 

U->U' = gJJgl, 

where (?l,r £ G g i 0 bai- In a model of the generalized hidden local symmetry with Gi oca i = SU(AT/)l x SU(AMr, the 
vector and axial-vector mesons are introduced as gauge bosons associated with the hidden local symmetry [ill lifll . 
Ggiobai x Giocai are spontaneously broken down to global symmetry H g i Q bai = SU(IV)v- The vector and axial- 
vector mesons acquire masses through the Higgs mechanism with spontaneous breaking of the GHLS. We define 
representatives £l,r,m hr the quotient space G g i 0 bai x Gi oca i/-Hgiobai such that U = Cl£m£r, and hidden local gauge 
bosons, Vl,r- These fields, £l,r.,m and Vl.r, are transformed under Ggiobai x Gi oca i as 

£l,r -*■ Cl,r = G l ,r (®)£l,r MSl.r > 

£m —> £m = G l (ie)£m (^)Gr (%), 

Vl,r m = g l,r(>)Vl,r m G[ R (ai) - idfj,G LtR (x)Gl R (x), (2.1) 

Gl,r( 2 i) G Giocai- The covariant derivative for £l,r,m are defined as 

l = - iV l m £l + 

D^r = <9 m £r - iV r m £r + 

m = m — + *£mRr/j, (2-2) 

where C M and 7are external gauge fields which are associated with the global symmetry Ggiobai- The covariantized 
1-forms, 

«l,r m = --D^l.rCI.r) ( 2 - 3 ) 


F ^’ K ^ = 9^143^ - ^Vl.Rm - i \ V \ l,r m , Vl.ri/], (2.4) 


and field strength, 
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are mapped under the gauge transformations as 

Sl,r^ - 

-*• Gl.r^Sl.RmGl.rW, 

Sm /1 - 

-»■ G h (x)aMij,Gl(x), 

f(l.R) _ 

HV 

- GL,n(x)Ffr R) Gl R (x) 


From these field and their transformations, the Lagrangian is constructed as seen in the following. 

As far as the lowest-order derivative terms are concerned, the gauge invariance and even-parity condition cast the 
Lagrangian into the following form in general; 


£( 2 ) = aF^trSy^ + 6F 2 trS^ 

+ cF 2 traM M + dF 2 tr(S A/ x + 5 m m ) 2 , 

( 2 . 6 ) 

where, 



S V/ x = ^ (CmSr^^m + “r.,i 

) ) “A /x = 2 (CmSr^^m “ “L/x) ■ 

(2.7) 

We can rewrite the Lagrangian by introducing a nii = Sa/x + 2m/x : 



£(2) = ( atr [ S v>] + - a MfJ ,) 2 } + ctr[aj^J + dtr[a 2 J) 

= F^(atr[Sy M ] + {b + c)tr[(a M/i - 7^— S^) 2 ] +{d + 


IT (1 

be . 

-)trfa 2 


7 Tfl 


])■ 


( 2 . 8 ) 


By renormalizing the pion field and appropriately rewriting the coefficients we obtain 

£(2) = -F 2 (tr[S 2 M ] + atr[S^J + /?tr[(a M/ x - 7^) 2 ]). ( 2 . 9 ) 

The first term is nothing but the non-linear sigma model Lagrangian. The second and third terms correspond to 
mass terms for vector and axial-vector mesons, n — A\ mixing is controlled by the parameter, 7 , which takes a value 
between 0 and 1. Unless the kinetic terms of the gauge bosons are incorporated, the Lagrangian is reduced to the 
non-linear sigma model as it should be in the low-energy limit. In fact, the equations of motion for V lt and A M are 
given by 


Vft = - (Cm m (o;v^ — otyiti + 7 a A/i)?M/j + avyj + a m m — 

9 - CXM 11 + 70!A/x)£m/x ~ % - + 7 a A/x) 


A,= 


( 2 . 10 ) 


where fields ay,A,M are defined as Sv.a.m with Vl,r = 0 . Substituting Ea. 12.101) into the Lagrangian Ea. 12.91) leaves 
only the first term: £( 2 ) = £ 2 tr[S 2 M ]. 

The kinetic term of gauge boson is given by 


£ ki n = -^tr[(i^)) 2 ] - ^tr[(fM) 2 ] 

= -^trl^W) 2 ] - ^tr[(F#>) 2 ], (2.11) 

where 

F™ = d^Vv - d v V^ - i[V„, K\ - i[A„, A„], 

KV = d»A v - d v A^ - i[V„,A v \ - i[A^V v \, (2.12) 

with = (Vr m + Ul/J /2 and = (Vr m - Vl m )/ 2 . 

Here we parameterize the nonlinear fields as 

£ R (ar) = 

^ h (x) = e ip(x) e ia{x) e-^ ix \ 

£ M {x) = e ip{x) e ipix) . 


(2.13) 




4 


By expanding the Lagrangian up to quadratic terms of these components, we obtain 

£quad =F 2 tr[(a(<9 /i er) 2 + (1 + /37 2 )(^7r) 2 + f3(d P p) 2 + aV 2 + /L4 2 

- 2 / 37 ^ 7 + 2/3(«9> - 7 ^tt)A m - 2ad p aV p 

- 2(3'yd p pAy - 2p 1 A p Ay + 2(1 + Ay + 200"^ - 2 aVyV 1 * + aV 2 + (1 + /3 7 V 2 ], (2.14) 

with V M = (JZy + Cy)/2 and = (JZy — Cy)/2. Note that the field 7r in this Lagrangian cannot be interpreted as the 
physical pion because 7r is mixed with p through the term — 2f3 r yd p ltd P p and components of 7r are partially absorbed 
by Ay,. To resolve the mixing and identify the physical pion, we shift and rescale the field as p —> {p/\J]3 + 77r)/F,r, 
which leads us to 


£quad =tr [(dya) 2 + (dy-rr) 2 + (dyp) 2 + M p p 2 + M Ai A\p + 2F p d p pAy - 2F a d p aVy 

- 2 F Al d»pAy - 2g/3 1 F 2 A p Ay + 2F^nA p + 2F a d p aVy - 2agF 2 V p Vy + F 2 V 2 + (F 2 + Fjj.4 2 ]. (2.15) 
Here, other fields are also rescaled as 


7r 


7r 

1 7 T 


a 


\faFn 


The gauge fields, decay constants, gauge couplings and mass of the gauge fields are expressed as 


(2.16) 


Pm = 9 Vy, 

Fa = 0.F 2 , 
Mp = ag 2 F 2 , 


Al M — 9 Ay, 

F 2 = JF 2 , Fl = f3^ 2 F 2 , 
M 2 Ai = !3g 2 F 2 . 


(2.17) 


III. WILSONIAN MATCHING FOR THE GHLS 


Since the GHLS model is an effective theory of QCD, the parameters of the GHLS model can and should be 
determined from QCD. In this section we apply the Wilsonian matching for the GHLS model [2 S|- Following 
ref. mm, we adjust the vector and axial-vector current correlators in the GHLS model to the correlators in QCD 
calculated by the operator product expansion (OPE). Here we assume the existence of an energy scale, A, above A\ 
meson mass where both the GHLS model and the OPE of QCD are applicable. 

The vector and axial-vector current correlators are defined by 


i J <Txe ipx ( 0 | TJ“(a;) J b ( 0 ) |0) = -5 ab P^n v (Q 2 ), 
* J d 4 xe* p * (0 \TJ^(x)JU0) | 0 > = -S ab P^H A (Q 2 ), 


(3.1) 


where Q 2 = — p 2 and P pv = q^p 2 — p p p u . These correlators are well described by the tree-level contributions around 
the matching scale, A |2J,[35|- In the GHLS model, the correlators up to 0(p 4 ) are given by 


n 


(GHLS) 


(Q 2 ) = 


F 2 (A) , Fi 1 (A)(l- 5 2 i (A)z 4 (A)/ 7 (A)) 


Q 2 

77 


M Al (A) + Q 2 


~2z 2 (A), 


tt(GHLS) //q2\ _ F 2 (A)(1-2 5 2 (A).~ 3 (A)) o 
v [Q) ~ M 2 (A) + Q 2 l(> ’ 


(3.2) 


where g 2 ( A) = 5 2 (A)/(1 + 2«(A)), g 2 Al (A) = ff 2 (A)/(l - 2«(A)), M 2 (A) = S 2 (A)F 2 (A) and A) = g Ai (A)F 2 (A). 


k. z i, z 2 , z 3 and Z 4 are coefficients in the following Lagrangian of 0(p 4 ): 


£(4 )z = Atr [F^ M F W ^l\ + 


Zi + z 2 


tr[(^) 2 + (PW 


Zi - Z2 


2g2 1 w SiV1 1 ^ IV 11V J \ [IIS J J 2 

+ ^ tr [( L ^)^(L) F + (R^)4f( R)F ] ' “ 3 "* 4 


tr [^ ) ^L^MCR^" (R)M!/ ^R.CM^L] 


(3.3) 
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On the other hand, the OPE in QCD leads to the correlators up to 0(1/Q 6 ) |3ftl37l as 

Q 2 , Tra a (G^G^) 14087r 3 a s ( qq ) 2 


n^ CD) (o 2 ) = 
n^ QCD) (Q 2 ) = 


l 

8^2 

1 

87T 2 


0 + t) 


In ■ 




Q 4 


81 


Q 6 


In 


Q 2 na s (G^G^) 8967r 3 a s (qq) 2 




Q 4 


81 


Q 6 


(3.4) 


where p, is the renormalization scale of QCD. 

Now let us make the Wilsonian matching of Eq. 13.211 with Eq. El at the scale A. We note that the correlators in 
Eg. !!!?"! explicitly depend on the renormalization scale p. To eliminate the explicit p dependence, we consider the 
difference, 11^ — Ily, and the first derivative of the correlators. Then the matching condition for the difference is 


(< GHLS) (Q 2 )-n(f HLS) (Q 2 ))| Q2 _ A2 = (n^ CD) (Q 2 ) - n^ QCD) (Q 2 )) 

The matching conditions of the first derivative for each n _4 y are given by 


Q 2 = A 2 


32-7rq 5 (qq) 2 
9 A 6 


-e 2 ^ <v HLS, (e 2 : 


Q 2 —A 2 dQ 2 


(3.5) 


(3.6) 


Q 2 = A 2 


where 




(GHLS) 


(Q 2 ) 


F 2 ( A) , A 2 F 2 i (A)(1-2 5 2 i (A)^ 4 (A)/ 7 (A)) 


-q 2 ^< hls V: 


Q 2 = A 2 


Q 2 = A 2 


A 2 


(A 2 + Af 2 i (A)) 2 


A 2 F 2 (A)(l-2 5 2 (A)z 3 (A)) 
(A 2 + A/ 2 (A)) 2 


(3.7) 


and 


- g *_Tn« CD) (<3 2 > 

- « 2 ^ni? CD> «5 2 ) 


Q 2 =A 2 


Q 2 =A 2 


1 

87T 2 

1 

8^2 


1 + 2 - 

7T 


1 +^ 

7r 


+ 


2tt a s (G^G^) , 14087r 3 a 5 (gg) 
3 A 4 

2™ s (G^G^ 


27 A 6 
8967r 3 q s (qg) 2 


A 4 


27 


A 6 


(3.8) 


When the chiral symmetry is restored, the vector current correlator coincides with the axial one due to (qq) —> 0: 


(n^ GHLS) (Q 2 ) - n[, GHLS) (Q 2 )) 


Q 2 = A 2 


0 . 


(3.9) 


The derivatives of the vector and axial-vector correlators of Eq. (13.81) become also identical with each other. We 
emphasize that their values are finite but nonzero because (a s / tt) (G fJ , v G ,lv ) may not vanish at the restoration point. 
This implies 


( F 2 i (A)(1-2 5 2 i (A)z 4 (A)/ 7 (A)) /n 

K(A) + - (i + m 2 Ai (A)/a 2 ) 2 -^ °’ 

F ff 2 (A)(l-2 ff 2 (A)^ 3 (A))^0. (3.10) 

With the parameters obtained by the Wilsonian matching, the low-energy parameters such as on-shell p meson mass 
are derived by the renormalization group. We discuss more detail of the restoration of chiral symmetry in seclvl 


IV. RENORMALIZATION GROUP EQUATIONS 

In the HLS model, quantum effects on couplings have been studied by the renormalization group [23l (iy| ■ It was 
pointed out that the quadratic divergence plays an important role in the chiral transition, especially in connection 
with the large Nf QCD |22]. Quantum corrections to the bare parameters suggest a violation of the vector meson 
dominance. Following this argument in the HLS model, we consider the renormalization group for the couplings in the 
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GHLS model to estimate quantum effects and find fixed points. To this end, we introduce dimensionless parameters 
for convenience; 


x(ii) = 


Nf M 2 
(4 n^F^Y 


G(n) = 


N. 


f „2 


(47t) : 


<? 2 (/4- 


(4.1) 


By calculating one-loop diagrams contributing to the bare couplings as shown in Appendix, we obtain renormalization 
group equations for the parameters as 


dx 


dp 


ay 2 /3y 2 


p— = x <2 — x \2 — a + 2y + 2 ay-----ay' 


P 


.2 _2 o„2.,2 , ft2 7 2 o 2 , /3 2 7 2 . 2 ,4 


+ -G a 2 - ay 2 - 2a 2 y 2 + 


P 


- Pi + 


da 

dp 


3 a 2 a 2 

~Y + 2p 2 


a 

, 2„,2 


+ a y 


2 2 , 2a2 7 2 , 3a 2 y 4 


/i— = £ { - - 2a + - 2ay - 3a y 


0 


+ ^Ga ( 1 - a 2 + ^ + 2ay 2 + 2a 2 y 2 - -/3y 2 - a 2 y 4 ) , 


dp 

dp 


a 

P 


p'2f- = x [2— 2- — 2/3 +a/3 — 2/3y 2 — 2a/3y 2 + ^ ^ + a/3y 4 


2„,2 


a 

2 a.,2 


3 2„.2 _ a2/37 4 


+ -G -3a + — + 4/3 - a /? - a/3y 2 + 2a 2 /?y 2 + 2/ry 2 - 


dy 2 2 

d-r- = 2xy 2 
dfi 

p d ^ = -G 2 

dfi 


p 

1 a 2 a ay 


n 2 / 1 a 2 2a a 2 2 a 2 y 2 

+ 3G 7 ( v - 1 -«^^2+^ + ^+«7 


44 1 

T~24 


2 , « 2 2 « o 2 O 2 2 , 2 ° 2 7 2 0(X .,2 , 2 P 2 1 2 


5 + a + 752-77^ 2a 7 - 2a 7 + 

P 2 P P 


- W + 


/3 

+ a z y' 


2^,4 


(4- 2 ) 


(4.3) 


If we take a limit of P —> 00 and y —> 0, the above RGEs for a; and a do not receive contributions from the Ai meson 
and are reduced to those in the HLS model. However, it is not the case for G, because the number of the gauge fields 
in the GHLS model is different from that in the HLS model and the loop contributions of the field p remain in the 
limit. This is due to the cancellation between the overall decay constants in the Lagrangian and denominators in the 
power series of the fields. In the next section we examine these equations for fixed points to discuss the properties on 
the chiral-transition point. 


V. RESTORATION PATTERN 

The chiral restoration is characterized by the condition that the on-shell pion decay constant vanishes, iG(0) —> 0. 
That is to say the non-vanishing of the dimensionless parameter x(p) defined in Ea. 14.11) in the low-energy limit: 

:r(0) y 2 0. (5.1) 

Inversely, x(0) = 0 means a broken phase because of F v { 0) y^= 0 and is excluded from the parameter sets associated 
with the restoration point. In addition to Eq. o. we analyze the restoration patterns of chiral symmetry using the 
vector and axial-vector current correlators. When the chiral symmetry is restored, the vector current correlator should 
be equal to the axial vector one, not only at the matching scale A as shown in Eq. cm but also at the low-energy 
region less than A. Accordingly, at the restoration point we require equality of these correlators; 

n v (Q 2 ) = n^Q 2 ), ( 5 . 2 ) 

for Q 2 < A 2 . Since the HLS model includes of only 7 r and p, the equality of correlators and Wilsonian matching 
uniquely bring about the VM scenario urn In the GHLS model, however, other possibilities exist. Actually 
considering the restraints on the decay constants of Eq. 13.11)11 . Eo. 15.211 gives Z\ (p) = Z 2 (p) and preserves the following 
three possibilities: 

(I) M 2 p (p) = M\ x {p), F\ x (p)( 1 - 2 g 2 Al (p)z 4 (p)) = F 2 (p){ 1 - 2g 2 p (p)z 3 (p)) and F 2 {p) = 0, (5.3) 

(II) M 2 {p) = 0, F 2 (p){ 1 - 2 g 2 {p)z 3 (p)) = F 2 (p) and F^p)^ - 2g 2 Ai (p)z 4 (p)/-f(p)) = 0, (5.4) 

(III) M p (p) = M Al (p) = 0 and F 2 (p) = F 2 (p)( 1 - 2g 2 p {g)z 3 (g)) - F\ x (p)(l - 2g 2 Ai (g)z 4 (g)/l(g)). (5.5) 
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In order to determine whether each of the above possibilities can be realized, we require that the restoration condition 
is satisfied for any energy below the matching scale. This demands that each of Eqs. 15. dll - 15.511 is renormalization- 
group invariant. 

Let us begin with the case (I). Eq. 15.51) . F 2 (p) = 0 implies that the power expansion of p 2 /(4 ttF 7I (A)) 2 breaks 
down. Therefore, the first possibility is inconsistent with the one-loop approximation adopted in this analysis and is 
not studied further. 

Next, we consider the case (II). From Eq. 15.411 we obtain 

g 2 = o, a = 1, /3 7 2 = 0, (5.6) 


for any p less than the matching scale. Since the renormalization group equations for these quantities are rendered as 


dG 

d-T~ = 0 , 
dfi 


da 1 


d(f3j 2 ) 

r 


= x 7 2 (-4 + - + 2 7 2 ), 


d(z ± - z 2 ) 


N f 1 

(4tt) 2 24/3 2 ’ 


(5.7) 


Ea, 15.(il) with z\{p) = z 2 (p) results in j3 —* oo and 7 = 0 for any p. As we take a low energy limit p —> 0, other 
parameters are also expected to be in a fixed point. Considering the above conditions together with the RGE for x 


dx 

p— = x{2 - x), 
dp 


(5.8) 


the non-trivial fixed point is found to be (x, a, /3 , 7 2 , G) = (2, l,oo,0, 0). This fixed point corresponds to the VM, in 
which M p —> = 0 and F a /Fn —> 1. We note that a combination of /3 7 2 remains to be zero as f3 —> oo and 7 —> 0 

from Eq. 15. (ill and hence Fa x = 0. We can say that A\ is decoupled to the system in this case and the GHLS model 
is reduced to the HLS model at the restoration point. 

Finally, we discuss the case (III). Note that the second condition of Ea. 15.51) is nothing but the first Weinberg sum 
rule which is saturated by the pole terms From Eq. 15.511 we obtain 

g 2 — 0 and a — /3 7 2 — 1 = 0, (5.9) 


for any p less than the matching scale. The RGEs for these quantities are given by 

dG d(a — /3 7 2 — 1) (1 — 7 2 )a 2 — 2a + 1 d{z\ — z 2 ) Nf (1 — 7 2 )a 2 — 2a + 1 

^ dp ’ ^ dp 2(1 — a) 2 ’ ^ dp (47t) 2 24(1 — a) 2 


(5.10) 


From these, we can see that a parameter set of a = 1/(1 — 7 ) and (3 = 1/((1 — 7 ) 7 ) makes Eci. 15.91) invariant under 
the renormalization group. Here we note that this parameter set satisfies the first and second Weinberg sum rules at 
0{p 2 ) even if we take a nonzero g |39). Under these conditions, remaining RGEs for x and 7 2 become 


P-j- = ^(2 - (! - 2 7 + 3j 2 )x), p^-j- 

dfi dfi 


= 2 7 2 (1 — 7 )(1 — 3 7 )x. 


(5.11) 


The parameters x and 7 are expected to be in a fixed point in the low energy limit of p —> 0. Except for the trivial 
fixed point of x = 0 with an arbitrary 7 , we find three fixed points from Eq. 15.11(1 . 


{ (1,1) (the standard type) 

(2,0) (the VM type) (5-12) 

(3,1/3) (the intermediate type). 

These fixed points are classified into the restoration types noted above according to the value of 7 , because the tt-Ai 
mixing are basically controlled by 7 as discussed in the next subsection. The parameters below the matching scale 
and the fixed points are summarized for each restoration pattern in Tabled 

Under the standard scenario, if a was in a fixed point at the matching scale, F a ( A) would diverge and be inconsistent 
with Eq. (l3~Tol because F a (A) must be finite. Thus we do not claim that the parameters are fixed on a fixed point at 
the matching scale. Under the VM scenario in the case of (II), the parameters are on a fixed point and 7 keeps taking 
zero for p lower than the matching scale. However, for the VM in the case of (III), 7 does not vanish for a non-zero 
p and is not on any fixed points. Consequently, in the case of (II) A\ is decoupled in the bare Lagrangian, while for 
the VM in case of (III) A\ remains to be dynamical at least near the matching scale. This difference is expected to 
reflect on the spectrum of A\ as the chiral symmetry is restored. 










0 < p < A 

fixed point (p —> 0) 


X 

a 

P 

7 

ip 

9 77 T 7 T 

(II) VM 

a = 1, /3 = oo, 7 = 0( /3 7 ' ! = 0 ) 

2 

i 

OO 

0 

0 

1/2 

(III) VM 
Standard 
Intermediate 

(a = 1/(1 — 7 ), 

\d = l/( 7 (l- 7 )) 

2 

1 

3 

i 

00 

3/2 

oo 

oo 

9/2 

0 

1 

1/3 

0 

tt /2 

arctan \[2 

1/2 

0 

1/3 


TABLE I: Parameters and their fixed points for the restoration patterns. The gauge coupling is zero, <7 = 0, for these restoration 
patterns. In the intermediate scenario, the fixed point is not realized in the limit of p —> 0 unless the parameters are adjusted 
at the matching scale. 


Although the parameters may not be located on the fixed point at the matching scale, the couplings must be on 
the critical line along which the renormalization-group flow enters into an infrared fixed point satisfying a;(0) / 0. 
The intermediate fixed point is on the critical line but an UV fixed point. It requires a fine adjustment to keep the 
point realized in IR region, but no reason exists to justify such an adjustment. The restoration conditions put severe 
constraints on the couplings, and without these conditions all of the fixed points become unstable. Unless Eq. O is 
satisfied, a fixed point of the vector realization defined as x = 0 with Eq. itO) will be achieved Eol. 

It turns out that the gauge coupling becomes zero, g = 0, at the all fixed points satisfying the restoration condition. 
This fact implies that p and Ai get lighter as the chiral symmetry is restored and eventually become massless at the 
restoration point. This result is consistent with the Brown-Rho scaling @. 


A. The mixing angle 

When mesons are classified by the light-front charge of the chiral group SU(Nf)L x SU(Nf)a in the broken phase, 
the spontaneous breaking of the chiral symmetry mixes the pseudo-scalar meson with the longitudinal part of the 
axial-vector meson while the scalar and vector mesons remain to be pure representations: 


|<7) = |(V 2 -l,l)©(l,iV 2 -l)}, 

K> = \(N f ,N* f ) © (N* f ,N f )) sin if) + \(N 2 f - 1 , 1 ) © ( 1 , 8 )) cos^, 

\p) = 

\Ai) = | (N f , N* f ) © {N},N f )) cos V -| (Nj- 1 , 1 ) © ( 1 , N] - 1 )) sin^. ( 5 . 13 ) 

The ex peri mental values of the decay constants suggest that the mixing angle, ip, takes ip ~ 7t/ 4 assuming two light 
flavors |23,[2a|- As the chiral symmetry is restored, the mixing is expected to disappear and physical states become 
pure representations because the light-front chiral charge commutes with the Hamiltonian j3(l ]. It is non-trivial to 
which representation the pion belongs. In the case that the pion belongs to ( Nf , N jf) © (iVjf, Nf) with ip = tt/2, the 
scalar meson becomes the chiral partner of the pion and we call it the standard scenario of the chiral restoration. In 
the case that the pion belongs to (V 2 — 1,1) © (1,-A 2 — 1) with ip = 0, the pion and the longitudinal mode of the p 
meson form a doublet, which is called the VM scenario ||. 

In this section, we determine the mixing angle at the restoration point for each restoration pattern classified in the 
above. We define the mixing angle in the GHLS model as the ratio of the pole residues for tv and Ai in the axial-vector 
current correlator. These residues correspond to tv and Ai decay constants squared and the mixing angle is expressed 
as 


tan 2 ip = 


^ 1 (A/i 1 )(l-gj 1 (M 2 i ) g 4 (M 2 i )/ 7 (M 2 i )) 

A 2 ( 0) 


/3(M 2 i ) 7 2 (M 2 i )(l-g 2 i (M 2 i ), 4 (M 2 i )/ 7 (M 2 i )). (5.14) 


At the restoration point, this expression becomes 


tan 2 ip = /3(0) 7 2 (0), 


(5.15) 


due to Ma x = 0 and g^i = 0. Therefore the mixing angle is controlled by the infrared behavior of the parameters 
near the restoration point. Corresponding to the three fixed points discussed in the previous section, the mixing angle 
at the restoration point is given by 

Standard :tan ip —> oo, VM : tanip —» 0, Intermediate : tan^ —> v2- 


(5.16) 
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As is expected, the mixing disappears at the standard and VM fixed points, but survives at the intermediate fixed 
point. The intermediate fixed point should be excluded if we require that the mixing of ir and A\ is resolved at the 
restoration point. 


B. Pion form factor and vector meson dominance 


The vector meson dominance (VMD) for the pion form factor is well established in the vacuum 0 . The pion form 
factor in the GHLS model is written up to 0(jp 2 ) as 


^ ± (p 2 ) = (l^^(l- 7 2 ))+^(l^ 7 2 ) 


Ml 


Ml — p 2 ' 


(5.17) 


The first term is the direct term, g 1T:w , and the second term is the p meson exchange term. The VMD implies 


= 1 - f (1 - 7 2 )« 0 . 

At the chiral restoration point, this coupling is expressed as 


(aO 


1 - 7(m) 
2 


(5.18) 


and its behavior at low energy is different according to the restoration patterns classified above. Under the VM 
scenario in the case of (II) g 77r7r takes a constant value of 1/2. In the case of (III), the coupling runs depending on 
the parameter 7 and shows (in)violation of the VMD in the infrared limit as 


0 

g-y 7T7T ' \ 1/2 

I 1/3 


(the standard type) 

(the VM type) as p —> 0. 

(the intermediate type) 


(5.19) 


In the standard scenario the VMD is satisfied. The violation of VMD appears in the VM and Intermediate fixed 
points. 

On the other hand, g pnn coupling is written by 


<W = fa(l-7 2 )- (5.20) 

The coupling g pnn becomes weaker as chiral symmetry is restored for all restoration patterns because the gauge 
coupling diminishes. This implies that p becomes stable for the decay into two pions near the restoration point and 
makes a sharp peak in the dilepton spectrum on the background from pions, universally. 

The difference among the restoration patterns may be reflected on the pion background through the violation of the 
VMD. The pion background in the dilepton spectrum is suppressed in the standard scenario, while the large violation 
of the VMD leads to an enhancement of the background in the VM scenario. Thus observations of the dilepton 
spectrum in the vector channel near the restoration point will ascertain the validity of the GHLS model adopted here 
as well as the true restoration pattern in the real world. 


VI. SUMMARY 

In this paper we have studied possible patterns of chiral restoration employing the GHLS model, which includes 
the axial-vector meson in addition to the pseudo-scalar and vector mesons. First, we obtained constraints on the 
parameters of the GHLS by matching the vector and axial-vector correlators of the GHLS model with those in the 
operator product expansion of QCD at a matching scale, A. Next, we derived renormalization group equations for 
the GHLS model by calculating one-loop corrections to the two-point functions. Then, we examined behavior of the 
parameters at low energy and searched for fixed points of them at the restoration point by the obtained renormalization 
group equations. As the restoration conditions, we required the equality of the vector and axial-vector correlators 
and the vanishing of the on-shell pion decay constant. 

We found four patterns of the running parameters at the restoration point which flow to three fixed points. The 
first fixed point correspond to the ‘standard scenario’, in which the vector and axial-vector mesons are degenerate. 




10 


The second one is the ‘vector manifestation scenario’, in which the longitudinal part of the vector meson and the 
pseudo-scalar meson are degenerate. In both these scenarios, the mixing of two representations associated with 7r and 
Ai mesons diminishes at the restoration point as expected. On the other hand, the third fixed-point is of a new type, 
which can be regarded as the ‘intermediate scenario’ in a sense that the mixing remains even at the restoration point. 
It is noted, however, that this fixed point is not an infrared one and will not be achieved unless the parameters are 
perfectly adjusted to the point at the matching scale. 

In each of the restoration patterns, both p and A\ are expected to become massless, since the gauge coupling 
vanishes at the restoration point. It is further expected that the vector meson becomes stable because the coupling 
constant, g P mr, also vanishes and the decay into pions is suppressed at the restoration point. Here we emphasize that 
these results are universal to all restoration patterns. These three scenarios, however, are different in the violation 
of the vector meson dominance. This difference will be reflected on the pion background in the spectrum of dilepton 
decay from the vector meson. Two restoration patterns running into the VM fixed point are also distinguished through 
the A\ spectrum near the restoration point owing to the different behavior of parameters at low energy. 

Now, a crucial question is whether these fixed points are physically realized or not. It is determined if one can 
reach the restoration points continuously as one approaches the critical temperature or density of chiral restoration. 
The application to the system of finite temperature and/or density is our future problem. 

We would like to mention here possible phenomenological consequences of the results of the present work. The 
dropping of the p meson mass implies the threshold enhancement in the p meson spectrum at some temperature 
and/or density because the p meson mass is expected to coincide with twice the pion mass, which is analogous to 
the enhancement of the sigma spectrum 0 in the standard scenario studied within the linear sigma model l l ‘A l43l| . 
(Here, explicit breaking of chiral symmetry is taken into consideration and the pion is supposed to have a small finite 
mass.) 

Another point to be mentioned is that the GHLS model does not include the scalar meson, which becomes the 
chiral partner of the pseudo-scalar meson in the standard scenario. To include the scalar meson in the analysis is 
therefore another important future problem. 

Acknowledgment 
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Note added 

After finishing up this work we became aware of similarly work by Harada and Sasaki p39| . Although their analysis 
takes a different parameterization and gauge fixing from ours, obtained renormalization group equations and results 
are consistent with our work. 


APPENDIX A: GAUGE FIXING 

We take a gauge fixing to calculate quantum corrections. We perform a R^ like BRS quantization. The BRS 
transformations are defined by 

<5b£v,a = *CV£v, a + *Ca£a,v, 

= i(Cy - Ca)£m — *£m(Cv + C A ), 

SbV p = d p Cy + i[Cv, V p ] + i[C A , A p ] = D p Cy, 

&bA p = d p C A + i[Cy, A p \ + i[C A , V/] = D P C A , 

S B Cy, A = *(CVCv,A + C A C A y), 

^B Cv,A = *Hv,A, 

d B H v , A = 0 , (A.l) 

where Cy tA and Cy.A are the Faddeev Popov ghost field, By_A is the Nakanishi-Lautrup (NL) field, and 

£v,a = -(£r ± Cl)- (A. 2 ) 

We choose the gauge fixing function as follows 

F = tr[Cv((?^V/ + ay Gy + —ay By)] + to\C A {d^ A p + a A G A + —a A B A )\, (A.3) 
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where 

a F 2 
3F 2 

Ga = -^WCa - £a + Cm-) - Cm-), 

Cm± = 2 (Cm ± Cm)- ( a - 4 ) 

We obtain the gauge-fixing and FP terms: 

£gf + £fp = — ifisF, 

Cqf = 2tr [Byid^Vfj, + a v G v + -ay-By)] + 2tr [B^d^Ap + oaGa + ^“a^a)], 

£ fp = 2itr[C v (d»DpC Y + a v G?V)] + 2itvC A (d^D^C A + a A G' A )], 
aF 2 

Gy = d B Gy = —^-(Gy Cv + CyGy + GaCa + CaGa), 

BF 2 

G' a = SbG a = —^ L (7(Gy Ca + CaGv + GaCv + Cv-Ga + GvCm- — Cm-Gv — GaCm+ — Cm+Ga) 

— (GvCm- — Cm-Gv — GaCm+ - Cm+Ga)). (A.5) 

By integrating B, £gf becomes 

Ggf = ——tr[(^V^ + a v (G v - ^-tr[G v ]) 2 ]-—tr [{d^Ap + a A (G A - ^-tr[G A ])) 2 ] 

tty iV f CCA tv f 

= tr[(0*%) 2 ] - — tr[(^A M ) 2 ] - 2tr[^^Gv] - 2tr[0M„G A ] 

ay a a 

- a v (tr[Gy] - ^-tr[G v ] 2 ) - a A (tr[G^] - ^-tr[G A ] 2 ), (A.6) 

where Nf is the number of flavor. 


APPENDIX B: RENORMALIZATION 

Counter terms from 0(p 4 ‘) Lagrangian are also necessary to cancel the divergences of one-loop diagrams in addition 
to counter terms from 0(p 2 ) Lagrangian. First, we define the renormalized parameters in the 0(p 2 ) Lagrangian 
denoted with the subscript r as: 

7T=Zy 2 7r r , a = Zl /2 <T T , p = Zp /2 p r + SnplTr, 

p = Zy 2 Pl , A, = Z^ 2 A lr , g 2 = Z g g 2 
F 2 — , F 2 = F 2 r + 5 F 2 , F 2 = F 2 t + S F 2 , F Ai = F\ iT + S F ^ , 

(B.l) 

where S^ p is necessary to resolve ir — p mixing in the loop level. The counter term is defined 5^ = 1 — Z^, where cf> = 7r, 
p. p , Aj, g. Z p = Z Al or 8z p = 8>z Al is satisfied by the parity invariance of the Lagrangian. We explicitly evaluate 
only the counter terms of 0(p 4 ) related to the vector and axial-vector current correlators. The relevant Lagrangian 
is given by 

F(a)z = ^2CmB*- R )^Cm] + 1+ ^ + (F^) 2 } + 1 9 2 tr[J r y ) Cj)CMCR^ r(R)M 1 / CRCMCL] 

+ fi ^ ii tr[CL + Z R rj£>& F W'»'] + li-iitr^yciCMB^^^ClrCL + ^CUL^^'CmCr]- 

(B.2) 


^■3,4 Zp Zg — £r3,4 T ^ 23 . 4 - 


We define renormalized parameters 

nZp = K r + S K , 


z l,2 — ~rl,2 + <>zi ,2 5 


(B.3) 
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t-ab TT/ii/ 


LL AA — 




+ 


+ 





FIG. 1: Contributions to II^ at one-loop. 


APPENDIX C: ONE-LOOP CORRECTIONS 

In this appendix, we evaluate some two point functions and their divergent part. We define A(m 2 ), B(p,mf,ml), 
B^(p,m 2 ,m%), B^(p,ml,m^), m \, m|) as follows: 


A(m 2 ) = J - 


d 4 k 


1 


(27t) 4 m 2 — k 2 ’ 
,4 7 


R , 2 2\ f d 4 fc _1_ 

!’ 2) J *( 27r )4 ( m 2 _ (jfc_p)2)( m 2 _ fc 2)’ 

D(1 , 2 2 n f d 4 &; & M 

B^(p, m 1 ,m 2 ) = / 4 . 2 —71- x 2u 2 —tar, 

J i{ 27 t)^ (mf — (k — p) z )(m2 — k ) 

nuv, 2 2 \ /' d 4 /c 

(P»m 1 ,m 2 ) - J . (2 ^ )4 _ ( k _ p ^( m 2 _ k 2y 

B^(v m 2 m 2 ) = f d4fc (2k-pY(2k-p) v 

° {jp ' m ^ m2 > J i( 2 TT)*(m 2 -(k-p) 2 ){m 2 ~k 2 )' 

The quadratic divergence plays an important role in the Wilsonian renormalization group equations. We identify the 
poles of D = 2 and D = 4 as the quadratic and logarithmic divergences after employing the dimensional regularization 
to take into account the quadratic divergence following Refs (23, denoted as: 


A = 


B = 


A 2 

(4tt) 2 ’ 


We evaluate the divergent part using A, B as follows 


(47t) : 


■In A 2 


(C.2) 


A(?n 2 )|div = A — m 2 B 1 
B{p,m\,ml) |div = B, 

B^(p, m 2 , m|)| d iv = y-B, 

B^{p,m\,ml )|div = - 9^(2 A- (m\ +m%)B) - P^y, 

P^(p,m 2 ,m|)| d iv = J^+PV^) 

= (-2A + (m 2 + m\)B ) - P^y + 

where = p 2 g^ - p»p v . 

We take the Landau gauge (ay = ola = 0) in the following calculations of one-loop corrections. 


(C.3) 
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1. A-A 

The contributions from each diagram shown in Fig. ^are 

(1) : f ') = -(1 + pY - a)N f 6 ab g^A( 0), 


( 2 ): 

(3) : 

(4) : 

(5) : 

( 6 ) : 


P 

P 


= - 9 2 F°;{a - pY) 2 N f S ab g^B(p, 0, M 2 ), 

= -g 2 F 2 ^ 2 N f 5 ab g^B(p,0,M 2 ), 

= -^p 2 'y 2 g 2 N f 6 ab g^B(p,0,Ml i ), 

= ^-'y 2 N f S ab B' J ‘ , '(p, 0 , 0 ), 

= - a N f S ab [(/3 7 2 - a) 2 B^(p, 0,0) - PY(PY - a)(B»(p, 0,0 )p v + p^B^p, 0,0)) 


(C.4) 

(C.5) 

(C. 6 ) 

(C.7) 

(C. 8 ) 


+ /3 2 7 4 pV5(p,0,0) 

With Ea. IIC.3ll the divergent parts are given by 

(1) :-7V /( 5“V"(l+/37 2 -«)A 

(2) : -N f 6 ab g^g 2 F 2 (a-P 7 2 ) 2 ^B, 

(3) : -N f 5 ab g^g 2 F 2 ^B, 

(4) : - Nf S ab g^g 2 F 2 ^^B, 


(C.9) 


(5) : ^ 7 2 N f S ab ^ (-2g^A - P^j + p^p v B^j , 

(6) : - a N f 5 ab ^(PY - a) 2 \ (-2g^A - + J(/3 7 2 + a)Vp^ 


(C.10) 


The sum of divergent parts is given by 


n AA\dW = N f 


g^ (A[ -(l + P'y 2 - a ) - ^y 2 - ^(/Ty 2 - a) 2 ') - ^9 2 F 2 B ( (a - PYf + PY + 


P^ v (p 2 1 


1 


P 


1 


- ^ ^ 7 2 + - (a - /W + JP^B ~Y + -(PY + ay 

12 \a a J 4 \a a 

The divergences can be canceled by the following counterterms: 

n^tldiv = -tfn + s n )<r - 2 S Z2 p^ + 8 AAP y\ 


(C.ll) 


(C.12) 


where we introduce the counterterm 5 AA to cancel the divergent part proportional to p^p". Such counterterm does 
not exist in the Lagrangian of any order. This divergence looks accidental. However this divergence is automatically 
canceled when one calculates an observable such as (T Au,(x)A I/ (0)). Therefore one may simply drop this divergence. 


2. V - V 

The contributions from each diagram shown in Fig. [3 are 

(1) : ( * ) = (1 + PY - a)N f S ab g^A( 0), 


(C.13) 
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s ab mx, = 





FIG. 2: Contributions to Ilyy at one-loop. 



The divergent parts are 


= |(«(1 - 7 2 ) - 2 ) 2 N f 6 ab B^(p, 0,0), 

= ±N f 6 ab B^(p, 0,0), 

= ^ N fS abB ^(p,0,0), 

= \~gN f S ab ((a - (3) 2 B^{ Pl 0,0) + pVP 2 B( P , 0,0) 
+ 2(o - 0)0(2p*B , '(p, 0,0) - 0,0))), 

= ~^^N f S ab g^B(p,0,M 2 Ai ), 

= - Fig- V(/3 - a) 2 N f 8 ab B{p, 0, M Ai ), 

= - aF 2 g 2 N f S ab g^B(p , 0, M 2 ). 


(1) :(l+/J7 2 -o)JV / 5°y“'4 

(3):±iV)* 06 (-2<TA 

(5) :^A^“ 6 ((a - /3) 2 (-2A<T - + /?Wb) , 

rp2 n 2 n 2 o 

(7) : - F 2 5 2 7 2 (/3 - a) 2 N f 6 ab g^B, 

( 8 ) : - aF 2 g 2 Nf5 ab g lll ''^B. 


(2) : —(a(l - 7 2 ) - 2 yN f 5 ab -2 g^A - P^ 


(C.14) 

(C.15) 

(C.16) 

(C.17) 

(C.18) 

(C.19) 

(C.20) 


(C.21) 
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FIG. 3: Contributions to 11^ at one-loop. 


The sum of divergent parts is given by 


n vvldiv = N f [<T (A ((1 + /3 7 2 - a) + =±(a( 1 - 7 2 ) - 2) 2 - i ^ - g(a - /3) 2 ) 
-^F 2 i?(^+ 5 2 7 2 (/?-a) 2 +a 5 2 
- ^B((a(l - 7 2 ) - 2) 2 + 1 + ^ + ^(a - /3) 2 ) + p^B 


(C.22) 


The divergent part is canceled by following counter terms: 


n vvldiv = - 2 S zl P^, 


(C.23) 


where we drop the divergent term proportional to . 


3. A-A! 


The contributions from each diagram shown in Fig. 0are 


(1) : (*') = —^Lg^i^a — (3 — P"/‘ 2 )N f S ab A(0), (C.24) 

''V' (i 

(2) : ■ } = Vl g ^N f 6 ab A( 0), (C.25) 

(3) : ( a ) = ^-^ v N f 6 ab A{ 0), (C.26) 

(4) : - (° \- = ^N f 6*(2(a-0)(jtf-a)B» v (p,O,O)-2(a + tf){a-fa)p> t p''B(p,O,O) 

V 7 T / za 4 

" ^ ^ 

+ 4(/3 7 2 - a)(a - /3 7 )p y 13 M (p, 0,0) - 4(a - (3){a + (3 r y 2 )p^'B v {p, 0,0) 

~2 ((/Ty 2 - a)(a - /3 7 ) - (a - /3)(a + /3 7 2 )) p^p^B^p, 0,0)), (C.27) 

(5) : ° V= = (2aB^(p, 0,0) - 2(a - (3'^p" B(p, 0,0) 

+ 4(a - f3 r ))p v B fJj {p, 0,0) - 4 apPB v [p, 0,0) + 2 fapl*p v B{p, 0,0)), (C.28) 

(6) : ~J^)= =aF^-a)g 2 N f 5 ab ^B{p,0,M 2 p )-^(B^{p,0,M 2 p )-B^(p,0,0))y (C.29) 


(7) : = -giF 2 (a - /3 7 2 )(/3 - a)g 2 N f S ab (g*"B{p , 0, M 2 ) - ^ (>>, 0, M 2 ) - B^(p, 0, 0))) . 

(C.30) 
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FIG. 4: Contributions to 11^7 at one-loop. 


The divergent parts are given by 

(1) :^ 5 ^(2 a-(3-^)N f 8 ab A, 

(2) : 21gi"'N f 5 ab A, 

(3) : 9 -l^N f 5 ab A, 

(4) :£jV>$“ 6 ((« - - a) (-2 g^A - P^|) - (a + /3 7 2 )(a - , 

(5) :£iV>$“ 6 i ( 2 « (-2^^ - P^|) - 2(o - /wys), 

(6) :fi-P 2 7(/3 - a)g 2 g^N f 5 ab B^ 

(7) : - <? 7 F 2 (a - /?7 2 )(/3 - a^g^NfS^B^ 

The sum of divergent parts is 


n AA! | div = N f91 


A 


P 1 


SV 77 -(2a ~ P~Pl ) + H-(a - /?)(/?7 - a) - 1 


a a 


+ B-Flg - 2 ((/? - a) - (a - /? 7 2 )(/3 - a))) 


/?)(/?7 2 - a) + 1) (a + /3 7 2 + l) 


12 Va 

The divergences are canceled out by the following counterterms: 


EC, 


_ //i/ 


i — /-». 

*4^1 Idiv * 


gF p F Al /H 


5p2 


+ 6z A1 -^+6z 9 )-2g5 Zi P^. 
Ai P P ) 


where we drop the divergent term proportional to p^p v . 


(C.31) 


(C.32) 


(C.33) 


4. V-p 

The contributions from each diagram shown in Fig. 0]are 


(1): 

/' 7r \ 

; 1 

= ~\g( 2C 2 - a 7 2 - 

(2): 

/'V'; 

= ^iV/<S“V^(0), 

(3): 

' 

; 1 

= ^gNf8 ab g^ 11 A{fi), 


(C.34) 

(C.35) 


(C.36) 
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'.= = o05(l - 7")(2 - o(l - 7 2 ))N f S ab B^(p, 0,0), 


(C.37) 


= - 5 lV /( ^^(p,0,0), 


(C.38) 


WL^±agN f 6 ab B^(p, 0,0), 


(C.39) 


: -^JV/J o6 (7(o - /3) 2 S Mi/ (p, o, 0) - /?(o - f3 1 )py i 'B(p, 0,0) 

- 2/3'y(a - f3)p^B v + 2(o - /3)(o - f3j)p u B 11 - (a - /?)(a - 2 /3'y)p lx p v B), 

■ FlgS\P - a) 2 JV /( y° 6 ^B(p, 0, M 2 J - ^ (>>, 0, M 2 J - 0,0 


(C.40) 


(C.41) 


y^(3-a)N f 6 ab (g^B(p,0,M 2 Ai ) - (>>, 0, M 2 J - P^(p, 0,0))) 


(C.42) 

(C.43) 


The divergent parts are given by 


(1) : - - 5 (2/3 7 2 - a 7 2 - a)N f 6 ab g^A, 

(2) ■\9^Nf8 ab g IM 'A, 

(3) :^7V /( 5“V^, 

(4) :±a 5 (l - 7 2 )(2 - o(l - 'y 2 ))NfS ab (-2g^A - P^|) , 

(5) +gN f 6 ab (-2<TA-P'“'|), 

(6) ■■^^agN f 6 ab (-2g^A-P^^] , 


79 !\T.xabf^v„ o \2 f O „LLV A 


(7) :—^N f 5 ab (j(a - /3) 2 ^-2^M - P^-J + /3(/3 7 - «)pVBj , 

(8) :gN f S ab g^ 2 (p-a) 2 ^F 2 g 2 B, 

(9) : - 9^(0 - a)N f S ab F 2 g 2 jB, 

( 10 ) : 0 . 

The sum of divergent parts is given by 

n vp = 9^/ g [a Q - /?7 2 + iay 2 + \ a + Jp + -J La ( 1 - 7 2 )(2 - a(l - y 2 )) 

-1 - {1 w^ a + k {a ~ /J)2 ) + \ fWb - °) 2 -§«*-«))) 


(C.44) 


- P^| (a(l - 7 2 )(2 - a(l - 7 2 )) + 1 + 


+ —j-ia-P) 2 ^ -p^p”B^ffi - a) . (C.45) 
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FIG. 5: Contributions to II 'i'i at one-loop. 


The divergence is canceled by the following counterterms: 

F a9 [Jn 


\\P- V 

U Vp 


div 2 y^ + ^ g +5 Zp )9^-2g6 Zi P^. 


where we drop the divergent term proportional to p^p”. 


(C.46) 


5. p- p 

The contributions from each diagram shown in Fig. 0are 


( 1 ): ; 


( 2 ): { 


V \ 


= -g 2 (a- /3)y 2 N f 5 ab g^A(0), 

= ^(a-mfS ab g^A( 0 ), 



= 0, 


= g 2 N f S ab g^ {D d 1)2 A(M 2 ), 
= g 2 N f S ab g^ {D d 1)2 A(M 2 Ai ), 


(6): =; i= =-aY(l- 7 fJV/^>,0,0) 


(7) : 

( 8 ) : 


= —(2(3-a) 2 g 2 N f 5 ab B^(p,0,0), 
= \g 2 B^{p, 0,0), 


(C.47) 

(C.48) 

(C.49) 

(C.50) 

(C.51) 

(C.52) 

(C.53) 

(C.54) 


/ p 


(9): =; 1= =— g 2 N f 6 ab ( 7 2 (a-f3) 2 B^(p,O,O) + (0 7 -a) 2 pVB(p,O,O) 

4 P 

+ 2y(/3 - a )(/?7 - a){p li B v (p , 0,0) + p^B^ip, 0,0) - p^p^Bip, 0,0))), 


(C.55) 
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( 10 ): 



( 11 ): 



= ~ a)Wg A N f 5 ab (g» v B(p, 0, M 2 m ) - (>>, 0, M 2 m ) - B^(p, 0,0)) ) , 

(C.56) 

= 0, (C.57) 



-y (/? -a) 2 g 4 N f S ab (g^B(p,0,M 2 Ai ) (b^(p, 0,M 2 Ai ) - B^(p, 0,0))) , 


(C.58) 


= b 2N > s "I 7 


d 4 fc 


i(27r) 4 {{k - p ) 2 - M 2 )(fc 2 - M 2 ) 
x {-2pV CT + W P + (p - 2fc)'V <T }{-2p p y <T ' + 2 p a 'g v P' + (p - 2fc) V V } 


fc p fc p / {k-p)a(k~p) c 

X (to- - —)(to- (t _ p) 2 


(C.59) 


= \g 2 N f 6 ab J 


d 4 k 


1 


i(27r) 4 (( k-p ) 2 - M Ai )(k 2 - M A J 
x {-2pV CT + 2/s w + (p - 2fc)' i ^ <T }{-2/(? iy<T ' + 2p a ' g up ' + (p - 2fc)V’ V } 


/ kpkp' . , 

x \9pp' ^2 )\9<y<y' ~ 


{k — p)a{k — p) c 
(,k — p ) 2 


= 0,0) - pVB(p, 0,0)), 


0 , 0 ) - 0 , 0 )), 


(C.60) 

(C.61) 

(C.62) 


where D is the space-time dimension. To keep the gauge invariance we take D —> 4 after the sum of diagram. 
The divergent parts are given by 

(1 ):-ff 2 (a-Ph 2 N f 6 ab g> u 'A, 

{2):^{a-p)N f 8 ab g pv A, 

(3) : 0, 

( 4 ) : g 2 N f S ab g p AD D 1)2 A(M 2 p ), 

(5) : g 2 N f S ab g p AD ^ A(M 2 Ai ), 

(6) : \a 2 g 2 (l - j 2 ) 2 N f S ab (-2g^A - P^j^j , 

(7) : ^(2/3 - a) 2 g 2 NfS ab (~2g^A - , 

(8) : ^9 2 Nf5 ab (-2g^ A - P^ 

(9) : (V(-2 7 2 (a - /3) 2 A + (/J 7 - a) 2 p 2 B) - P pv + ^ _ a) 2 ^ ^ ; 

(10) : -F 2 (P-a) 2 'y 2 g 2 g 2 N f S ab g^B, 

( 11 ) : 0 , 
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(12 ):-|(^-a)VsW 1 '^, 


1 


,14 


(13) : -zgNfS P^-B - 2{D - 1 )<f v A[Mi) - 


2 , P 2 9^ 


B 


(14) : \g 2 N f 5 ab (VyS - 2 (D - l)g^A M ^ - 

(15) : -^-N f 6 ab (^-g^{2A + p 2 B) + , 

(16) : -^N f 6 ab (-g^(2A + p 2 B) + P^p^j ■ 

The divergent part including the internal line of p is given by 

(D — l) 2 


P 2 g^ 


B 


(4) + (13) + (15) = g 2 N f S ab gT 


D 


A(M 2 ) - (D - 1 )A(M 2 ) -^jB+ 1(2 A + p 2 B) 


+ P ,1V ( -B~- — 

3 4 3 


= g 2 N f S ab g* v 


1 


1 


(1 - -D)A(M 2 ) - -A(M 2 ) + -A) + P^—B 


,13 


6 


= g 2 N f S ab g' lv -M 2 B + P^ — B 
y : y 4 p g 


where we take D — > 4 and use 


lim —(1 - -D)A(M 2 ) = -M 2 B. 
d^aD 2 K pJ 4 p 


Similarly the divergent part including the internal line of Ai is 

(5) + (14) + (16) = g 2 N f S ab 
The sum of divergent part is given by 

n wldiv = 9 2 N f 5 ab 


g ,JV -Ml B + P^—B 
y A Al 6 


(C.63) 


(C.64) 


(C.65) 


(C.66) 


-(« - 0 T 2 + -J (o - 0 - \« 2 0 - - 7 2 ) 2 - ^( 2/3 - a ) 2 - \ + ^ ^ 


- F' 2 g 2 -B ( (/3 - a) 2 -l(/3 7 2 + 1) - (a + /3) 


+ P IJV B ( y + y (y« 2 (l - 7 2 ) 2 - y (2/3 - a) 2 - 1 - 


2 y 2 (a —/3) 2 6 


0 




The divergent part is canceled out by the following counterterms: 


n pp Idiv = -( M P 2 (^ S + S Zp ) + g 2 S F 2 )g^ + (5 Zp + 2S K )P^. 


(C.67) 

(C.68) 


6. Ti - Ai 

The contributions from each diagram shown in Fig. [U]are 


(1): 

/ 7T \ 

i j 

= g 2 (a-/3) 7 2 N f S ab g^A(0) 

(2): 

A'p'\ 

i * 

= pa-/3)N f S ab g^A(0), 


(C.69) 


(C.70) 
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n AiA, - 


+ 



+ 




FIG. 6: Contributions to II^ Ai at one-loop. 


(3): 



= 0 , 


(4): 



(5) : 

( 6 ) : 



(7): 



( 8 ): 



= g 2 N f S ab g^ {D D 1)2 A(M 2 p ), 

= g 2 N f d ab g^ v (D D 1)2 

= ^ N f S ab ((a - P) 2 B^(p, 0,0) + (0 7 - afp p p v B(p, 0,0) 

+ 2(0 - a)(/3j - a)(jfB v (p, 0,0) + p v B p {p, 0,0) - 0,0))), 

= £pN f 6 ab (a 2 B^(p, 0,0) + (/3 7 - a) 2 p^B(p, 0,0) 

+ 2 a(a - (3' 1 )(j?B v {p, 0,0) + p v B p {p, 0,0) - //5(p, 0,0)), 

= 0 , 



( 12 ): 



= -F 2 {P - a)Wg A N f 5 ab (g» v B(jp, 0, M 2 ) - ± ^(p, 0, M 2 ) - B^{p, 0,0)) ) 
= ~ (0- a ) 2 g A Nf5 ab (g^B(p,0,M 2 p ) - ^ (>>,0,M 2 ) - B^(p, 0,0))) , 


= g 2 N f S ab J 


d 4 fc 


i(2n)* ((k - p)* - M*)(k 2 - M\J 

x {-2 p p g p(T + 2 p a g pf> + (p - 2fc) V <T }{-2p'V' T ' + 2p <T V P ' + (p - 2fc) V V } 

„ {k - p)a{k - p)„'^ 

x ti/pp' , 2 )\y<Tcr' /, n 2 /’ 


fc 2 


{k-pf 


= -j-N f 5 ab {B p "(p, 0,0) - 0,0)). 


The divergent parts are given by 

(l):g 2 (a-Ph 2 N f 6 ab g^A, 
(2 ):^(a-p)N f S ab g^A, 
(3) : 0, 


(C.71) 

(C.72) 

(C.73) 

(C.74) 

(C.75) 

(C.76) 

, (C.77) 
(C.78) 

(C.79) 

(C.80) 
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FIG. 7: Contributions to 11,™ at one-loop. 


(4) : g 2 N f 6 ab g^ {D d 1] \ a - M 2 p B), 

(5) : g 2 N f 6 ab g^ {D d 1] \ a - M 2 Ai B), 

(6) : ^-N f S ab ^-2. \g^(a - f3) 2 A - P^{a - (3) 2 B + (/3 7 - afp^B^j , 

(7) : w Nf5ab ( ~ 2a2Ag ^ - p ^i B + (/?7 -°)w). 

(8) : 0, 

(9) : —Nf5 ab F 2 (f3 - afPg^^-^B, 

(10) : -N f 5 ab ^ ((3 - afg^^^B, 

(11) : N f 6 ab g 2 Ip^—B - 2 g^A + 3 g^(M 2 + M 2 Al )B - ^—^—B \ , 

(12) : -y N f S ab (^-g^(2A + p 2 B) + P^^pj . (C.81) 

(4) + (5) + (11) + (12 ) = g 2 N f 5 ab g^(M 2 + M 2 Ai )B + . (C.82) 

The sum of divergent parts is given by 

| div = S 2 ^ [<T ((« - /?)7 2 + £(« " 0 " ~ 

- \ F W B ((P - a) V + ^(/? - a) 2 - (a + /3))^ 

+ P^B f ■y - ^-(a - /?) 2 + 3^(/3 7 - a) 2 + T« 2 + 3^(/3 7 - a) 2 )] . (C.83) 

\ 3 12 a a ap ap ) 

The divergent part is canceled out by the following counterterms: 

IIZa, | div = -(M 2 Ai (S Zg + S ZAi ) + gH n )g^ + {5 Zp 2 S K )P^. 


(C.84) 
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7. 7T — 7T 

The contributions from each diagram shown in Fig. Qare 
'' n \ 4 — 3a(l — 7 2 ) 2 p 2 


( 1 ): 

( 2 ): 

(3) : 

(4) : 

(5) : 

( 6 ): - 

(7) : 

( 8 ) : 

(9) : -J 

( 10 ) : -J 

( 11 ) : 


12 


= 


F 2 


N f S ab A( 0), 




= %^ ab m, 


= —g 2 (a - f3h 2 N f 6 ab (l - D)A(M 2 ), 

= g 2 (a-P)'y 2 6 ab N f (l-D)A(M 2 Ai ), 
a(l — 7 2 ) 2 p 2 


FI 


N f 6 ab (A(0)+p 2 B(p,0,0)), 


2 2 
7 P 


= ^ j^N f 6 ab ((-4(3a - 2 (3 2 + a 2 )+(0) + 03- a) 2 p 2 B(p , 0,0)) 


p 

p 


= -<?V(1 - 7 2 )^ /( 5“ b ( p 2 B(p, 0, M 2 ) - 
-1 


2\ PV 


P ' M 2 


= -jg 2 (a - P) 2 ^ 2 N f S ab ( p 2 S(p, 0, M 2 ) - 


2 \ PV 


P ' M 2 


-7 2 2, 


g 2 {p-a) 2 N f 5 ab p 2 B(p,0,M 2 J- 


f 2,_P_V 

Ml 


(C.85) 
(C.86) 
(C.87) 
(C.88) 
(C.89) 
(C.90) 
(C.91) 

(^(p,0,M 2 )-B^(p,0,0))) , (C.92) 

(b mv (p,0,M2)-B mi/ (p,0,0))) , (C.93) 

(B^(p,0,Ml 1 )-B #1I/ ( P) 0 ) 0))^ , (C.94) 



- = F 2 g 4 (P - a) 2 1 2 Nf 5 ab [{D - l)2?(p, M 2 , M 2 J - ^ (B(p, M Ai ,M 2 ) - £(p, 0, M 2 )) 

+ m * ] ~ hilv{p ’°’ M * )} )■ (c - 95) 

The divergent parts are give by 


0) : * 

(2, :^AN t ^ A , 

a F 2 

2 i 2 

,-ox 7 1 P M cab a 

f? ’ ’ 

(4) :g 2 (a - (3) 1 2 N f 5 ab {D - 1)(A - M 2 S), 

(5) : - g 2 (a - p)j 2 N f 6 ab (D - 1)(A - M\B\ 

(s) a A^rA N ,^ (A+V 4 B) , 

(7) - 2/? 2 ^ 4o/?)24 + (/? - a) 2 p 2 B), 

(8) :-p 2 g 2 a 2 (l-^) 2 N f S ab ^B, 
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(9) '■—^-g 2 Nf8 ab (a — /?) 2 7 2 |.B, 

(10 ):ZPhl g 2 Nf s a \p- a )^ Bl 

(11) :F 2 g\f3 - a)WN f 5 ab {D - 1 )B. 


The sum of divergent parts is given by 

" 2 r / 4- 3a(l - y 2 ) 2 


II I — —Nf — 

ll7r7r ldiv — f p2 


A - 


12 


Pr 

a 


7 2 1 a( 1 — 7 2 ) 2 7 : 


ji“ + 


1 1. 


(3 a' 


+ 4^'^ + “ 2) 


( « 2 (! ^ 7 2 ) 2 + 7 2 (« - + -) ) + V s ( «(1 - 7 2 ) 2 + ^ ~ a) 


af3 


The divergences are canceled out by the following counterterms: 

H-7T7T I c jJ v P T V ^47r7r); 

where S^ nw is the counter term from the 0(p 4 ) Lagrangian. 


(C.96) 


(C.97) 


(C.98) 


8. 7T — p 

The contributions from each diagram shown in Fig. IHlare 


(1) : 

/' 7r \ 

1 


(C.99) 

(2): 

(T \ 

1 | 

= 2^aFf fS ‘ tA( ° ) ' 

(C.100) 

(3): 

1 l , 

2v^3 6/1 F 2 f 1 

(C.101) 

(4): 

O 

= ^=2g 2 (a- /3)3N f 5 ab A(M 2 ), 

(C.102) 

(5): 

O 

= ^=2g 2 ((3~a)3N f 6 ab A(M 2 Ai ), 

(C.103) 

(6): 

. " " N 

/ <T \ 

-1 r — 

\ K y 

= - l 2 )^N f 5 ab ((a + 2/3)71(0) + (a + /3)p 2 l?(p, 0,0)), 

(C.104) 

(7): 

/'p~\ 

r— 

\ ° J 

= - 2/3)71(0) + (a - /3> 2 l?(p, 0,0)), 

(C.105) 
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( 8 ) : 


(9): -J 


( 10 ) : 


( 11 ) : -- 


^=2/a(l - 7 2 )(a - (3)Nf^ ab [p 2 B(p, 0, M 2 ) - ^ 


_]_ 2g 2(2l_o)_ i a_ f3 )NfS ab (> B(p,0,MD - 


7 o„2 


2g 2 (f3-a)N f 6 ab (p 2 B(p,0,M 2 p )- 


r2\ P # V' 


P’ M2 


) - i|r (b^p,o,m 2 ) -b^(p,0,0))) , 

(C.106) 

I - (b^(p, 0, M 2 ) - B pu (p, 0,0)) ) , 

(C.107) 

(^(p,0,M2)-B^(p,0,0))^ , (C.108) 


(C.107) 



fg 4 Fl(a - P) 2 N f S ab (^3B(p, M 2 p , M 2 Al ) - (S(p, M 2 i; M 2 ) - £(p, 0, M 2 )) 


+ jyj-2 ^p-2 (B[iv{p, M Al , M p ) B pv (p, 0, M p 


(C.109) 


The divergent parts are given by 


7 2/3 — 3a 


^( 1 - 7 2 )| zNjS^A, 


^ 2 7?-ai Nf6abA ' 

(3) -w^ri NfSabA 

( 4 ) :-^j=2g 2 {a - f3)W f 5 ab (A - M 2 p B), 

(5) ^2 5 2 (/3 - a)3N f 5 ab (A - M Ai B), 

(6) :^^(1 - 7 2 )||iV/<5 o6 ((a + 2/J)A + (a + /3)p 2 i?), 

(7) : " W)A + (a -^ 2i3) ’ 

(8) :^=<? 2 a(l - 7 2 )(a - P)N f 6 ab p 2 ^B, 

( 9 ) :g 2 ^j^-(a-P)N f 6 ab p 2 ^B, 

(10) :0=((3-a)N f 6 ab p 2 ^B, 

(11) :^=2g 4 Fl(a- p) 2 3N f S ab B. 


(C.110) 


The sum of divergent parts is given by 

TT I _ 1 P Ar tab I" A f t 1 a 1 


2 vW 


P , -1,1, 1 


^N f S ab A { (-/3 - -a)(l - 7 2 ) + 7 + —(^ ^ o°) - (! - 7 2 )^(« + 2/3) + — (a - 2/3) 


a /3 3 2 


+ 2Fl-g 2 B —a(l- 7 2 )(/3-a) + 


(2/3 - a) 


(a - /3) + (/3 - a) 


+^-i3(-(l - 7 2 )(a +/3) + —(a -/3)) . 

The divergences are canceled out by the following counterterms: 

n^pldiv — V ( — fi-irp T P $4-71 -p)) 


(C.lll) 


(C.112) 


where S^p is the counterterm from 0(p 4 ). 
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APPENDIX D: RENORMALIZATION GROUP EQUATIONS 

The renormalization group equation for the parameter A, = F%,a,/3 ,..., at one loop is given by 


M dp dA 


(D.l) 


where e)A, is counter term for A;. From Eas. llC.12L HC.23L (IC.33L (IC.46L (IC.68L (IC.84L (1C.9811 and (IC.112L we 
obtain the renormalization group equations as follows: 


dF 2 

2N f 



(47t) 2 




+ 

dF 2 

2Nf 



(4n) 2 

d" 



+ 

dFl 

2N f 

2 

» dp - 

(4t r) 2 7 

< 

2/V f 


M dp 

(47t) 2 

d 



+ 

dg 2 

9 4 ^y 


y (4tt, 

dz\ 

/Vy 

1 

^ dp 

(4tt) 2 24 

dz 2 

Ay 

1 

^ dp 

(4tt) 2 12 

dz 3 

Ay 

1 

^ dp 

(4tt) 2 12 

d^4 

Ay 1 


(47t) 2 6 7 

dn 

d~r = 

Ay g 2 


2 


1 a 2 

4 + T 

p2 „2 


a 7 2 

/3 7 2 

a 7 4 \ 

2/3 

2a 

2 J 


„ 2_,2 


2a 2 7 2 + 

ft 7 

a 

-/?7 2 + 


aw , 2 4 

-hay 


a 

w- 

a 2 


a 2 7 2 a 2 y 2 

~Y~ + ~ 2 /T 

2 2 | /o2 _,2 


^T 2 aV 

2 + 4 


'2 / i a a a') 2 /?V 

{- 1 ~ a + 2a + a + — + wr 

+ \ f *9 2 (^a + 2a 2 - - 2 a/3 - a 2 y 2 + /? 2 q 2 + 2/3 

a ay 2 a 2 "/ 2 
2p- — + 2a 
a 2 

H—-—|- 4/3 

P 

22 1 


_ . _ 2 a 2 2a 2 2 o 2a 2 q 2 2 2/3“q 2 2 4 

y~ 481 5+ “ + ^-y - 2 “ 7 - 2 “ 7 + — - Vl + — + “ 7 


ft 


5 — 4ft + ft + —- — 2ft 7 + 
P z 

9fl 2 , /37 2 
a - 2/3q H- 


2a 2 y 2 


+ 2/3q 2 + a 2 y 4 


a 

" 2 2 a 


1 + 2a — a" — —yr H———h 2a7 + 2a 7 — 

a a 
aw 


2a 2 ') 2 


— 2/37 2 — a 2 7 4 


dp. 


2 a . 6a 2, 0 22 

-5- a - 7^ + yr + 607 + 2a 7 - 

a a 


2_,2 


2a 2 7 

“3“ 


- 6 / 3 7 2 + 


2 / 3 2 7 


2^,2 


2 4 
ft 7 


We define the parameter to calculate the RGEs for the dimensionless parameter for convenience: 


x{n) 

G(p) 


Nf P 

(4 WF 2 ' 

Nf g 2 . 

( 4 t r ) 2 ' 7 


(D.2) 


(D.3) 
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The beta functions for the dimensionless parameters x, a, /?, 7 2 , G are given by 


Px 

Pa 

P/3 

Pry 2 

Pg 




l 


dF? 


dFl 


F 2 ^ dg M dg^' 


1 


O 




PF% dn 
Nf dg 2 
(47 t) 2 ^ dg 


dF 2 

ar p 2 \ 

dfi 


We obtain the RGEs in Eq. (TOli by combining Eqs. cm and ID. 411 . 


(D.4) 


[1] T. Hatsuda and T. Kunihiro, Phys. Rept. 247, 221 (1994). 

[2] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000). 

[3] G. E. Brown and M. Rho, Phys. Rept. 363 , 85 (2002). 

[4] G. Agakichiev et al. [CERES Collaboration], Phys. Rev. Lett. 75, 1272 (1995). 

[5] F. Bonutti et al. [CHAOS Collaboration!, Phys. Rev. Lett. 77 , 603 (1996); F. Bonutti et al. [CHAOS collaboration!, Nucl. 
Phys. A 677 , 213 (2000). 

[6] K. Ozawa et al. [E325 Collaboration], Phys. Rev. Lett. 86, 5019 (2001); M. Naruki et al, arXiv:nucl-ex/0504016 

[7] E. Scomparin, Plenary talk at Quark Matter 2005. 

[8] M. Harada and K. Yamawaki, Phys. Rev. Lett. 86, 757 (2001) 

[9] G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 (1991). 

[10] M. Bando, T. Kugo, S. Uehara, K. Yamawaki and T. Yanagida, Phys. Rev. Lett. 54, 1215 (1985). 

[11] M. Bando, T. Kugo and K. Yamawaki, Nucl. Phys. B 259, 493 (1985). 

[12] M. Bando, T. Kugo and K. Yamawaki, Prog. Theor. Phys. 73, 1541 (1985). 

[13] T. Fujiwara, T. Kugo, H. Terao, S. Uehara and K. Yamawaki, Prog. Theor. Phys. 73, 926 (1985). 

[14] M. Bando, T. Fujikawa, and K. Yamawaki, Prog. Theor. Phys. 79, 1140 (1988) 

[15] C. N. Yang and R. L. Mills, Phys. Rev. 96, 191 (1954). 

[16] J. J. Sakurai, Currents and Mesons, Chicago University Press, Chicago, (1969). 

[17] J. S. Schwinger, Phys. Lett. B 24, 473 (1967). 

[18] J. Wess and B. Zumino, Phys. Rev. 163, 1727 (1967). 

[19] S. Weinberg, Phys. Rev. 166 , 1568 (1968). 

[20] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 

[21] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321, 311 (1989). 

[22] M. Tanabashi, Phys. Lett. B 316 , 534 (1993). 

[23] M. Harada and K. Yamawaki, Phys. Rev. Lett. 83, 3374 (1999). 

[24] M. Harada and K. Yamawaki, Phys. Rev. D64, 014023 (2001) 

[25] M. Harada and C. Sasaki, Phys. Lett. B 537, 280 (2002). 

[26] M. Harada, Y. Kim and M. Rho, Phys. Rev. D 66, 016003 (2002). 

[27] F. J. Gilman and H. Harai, Phys. Rev. 165, 1803 (1968) 

[28] S. Weinberg, Phys. Rev. 177 , 2604 (1969) 

[29] S. Weinberg,Phys. Rev. Lett. 65, 1177 (1990). 

[30] K. Yamawaki, arXiv:hep-th/9802037 

[31] M. Bando, T. Kugo and K. Yamawaki, Phys. Rept. 164, 217 (1988). 

[32] STAR Collaboration, Phys. Rev. Lett. 92, 092301 (2004) 

[33] G. E. Brown, C. H. Lee and M. Rho, arXiv:nucl-th/0507073 

[34] M. Harada and K. Yamawaki, Phys. Lett. B 297, 151 (1992). 

[35] M. Harada and K. Yamawaki, Phys. Rept. 381, 1 (2003). 

[36] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[37] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 448 (1979). 

[38] S. Weinberg, Phys. Rev. 177 , 2604 (1969). 

[39] M. Harada and C. Sasaki, arXiv:hep-ph/0511312 

[40] H. Georgi, Phys. Rev. Lett. 63, 1917 (1989). 

[41] T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55, 158 (1985). 

[42] S. Chiku and T. Hatsuda, Phys. Rev. D 58, 076001 (1998). 









28 


[43] Y. Hidaka, O. Morimatsu, T. Nishikawa and M. Ohtani, Phys. Rev. D 68, 111901 (2003). 



